Vibrational modes identify soft spots in a sheared disordered packing 
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We analyze low-frequency vibrational modes in a two-dimensional, zero-temperature, quasistat- 
ically sheared model glass to identify a population of structural "soft spots" where particle rear- 
rangements are initiated. The population of spots evolves slowly compared to the interval between 
particle rearrangements and the soft spots are structurally different from the rest of the system. 
Our results suggest that disordered solids flow via localized rearrangements that tend to occur at 
soft spots, which are analogous to dislocations in crystalline solids. 



Like liquids, solids can flow under applied shear 
stresses. Crystalline solids flow via rare rearrangements 
controlled by a population of lattice defects, namely dis- 
locations (jj. In disordered solids, rearrangements tend 
to be localized but there is no obvious way to iden- 
tify defects that might control them [IjQ. Can these rare 
localized rearrangements occur anywhere, as in a liquid, 
or do glasses possess a population of "soft spots," analo- 
gous to dislocations in crystalline solids, which are struc- 
turally distinct and susceptible to rearrangement? Al- 
though useful continuum models assume the latter [s, '?] , 
such a population of spots has never been identified from 
structural information. 

In order to search for a population of soft spots, we 
must start with a solidlike description of the glass. We 
begin with harmonic theory, in which the linear response 
to an applied stress is completely characterized by the 
normal modes of vibration. This approximation breaks 
down before solids begin to flow, so one would not expect 
the linear response to yield much insight into particle re- 
arrangements. However, recent evidence suggests that 
low-frequency vibrational modes, which are generically 
more prevalent in disordered solids than in crystalline 
ones |8l-[l0|. can be quasilocalized. Such modes have un- 
usually low energy barriers to rearrangement s 11 ill, and 
therefore are correlated with rearrangements |12h15| . 

In this paper, we use low-frequency modes to identify 
a population of soft spots in a model glass. We find that 
rearrangements begin at soft spots, the population of soft 
spots evolves slowly compared to the time between rear- 
rangements, and that there are structural differences be- 
tween soft spots and the rest of the system. We therefore 
conclude that soft spots are good candidates for elemen- 
tary defects that control the flow of disordered solids. 

We study a 50:50 binary mixture of soft discs of diame- 
ter ratio 1.4 in two dimensions, interacting via a Hertzian 
potential V = e(l — r/R)^/"^, where r is the distance be- 
tween the centers of two particles and R is the sum of 
their radii. Results presented here are for jammed pack- 
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ings with a packing fraction (j) — 0.95, which is much 
higher than the jamming transition at (fic — 0.84. We 
have also identified soft spots at values of closer to the 
transition jl6i] . Lengths and frequencies are in units of 
the small particle diameter and the interaction energy e. 
We employ Lees-Edwards boundary conditions to shear 
the system with a strain step of 10~^. After each strain 
step we relax the structure to its minimum energy to 
shear the system athermally and quasistatically. 

With increasing applied strain, the shear stress in- 
creases with a slope given by the shear modulus, then 
drops abruptly when there is a rearrangement. The 
strain step size is reduced to 2 x 10~^ before each re- 
arrangement. Between rearrangements, the dynamical 
matrix M is calculated at small strain intervals to ob- 
tain its eigenvalues (corresponding to the square of the 
frequency) and eigenvectors (the vibrational modes) [iTj . 

In this limit of zero temperature and strain rate, a 
rearrangement occurs when one vibrational mode (the 
critical mode) reaches zero frequency. At that critical 
strain, 7c, the packing becomes unstable and the coor- 
dination of particles in the packing changes. The initial 
rearrangement can trigger an avalanche of additional par- 
ticle motions [l^, so that the net displacements of the 
particles may be very different from the critical mode and 
may involve contributions from a number of modes [Tsf . 
However, in solids at finite temperatures and strain rates, 
fiuctuations can interrupt or extend avalanches. There- 
fore, for the remainder of this paper we focus not on the 
avalanche but on the reproducible initial particle rear- 
rangement, described by the critical mode [l8| . 

As the system is strained, the packing becomes less sta- 
ble and the mode frequencies tend to shift downwards. 
At a given strain, one might expect the lowest frequency 
mode to be the one whose frequency vanishes at the 
next rearrangement. Fig. [T] shows this is not generally 
true - the mode most similar to the critical mode lies 
at the lowest frequency only for a small range of strains 
immediately preceding the particle rearrangement [l5|. 
Note that most excitations at low frequencies are weakly- 
scattered sound waves with a strong plane-wave char- 
acter. These excitations coexist in the same frequency 
range as the quasilocalized excitations [ll|; as a result. 
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FIG. 1: (color online) The lowest ten normal mode frequen- 
cies as a function of applied strain. There are two critical 
strains at which a mode frequency approaches zero and par- 
ticles rearrange, at 7c — 0.014 and 7^ — 0.0183. The color 
of each point indicates the overlap of that mode with the 
critical mode at 7c — 0.0183. The lowest energy mode does 
not resemble the critical mode until just before the particle 
rearrangement . 
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FIG. 2: (color online) Soft spots in a system calculated at 
3.2 X 10"'^ units of strain before a particle rearrangement, (a) 
Regions of large displacement in theA^m = 30 lowest frequency 
modes. Bold circles highlight the Np = 20 particles with the 
largest polarization vectors, and different colors correspond 
to different modes, (b) Soft spots generated by clustering the 
particles highlighted in (a). Red arrows indicate the displace- 
ment of each particle during the next rearrangement. 



the normal modes exhibit characteristics of each. We 
therefore look at the entire population of low-frequency 
modes to extract soft spots, as follows. 

For a granular packing of = 2500 particles, we first 
identify the Nm lowest frequency modes in the spectrum 
of the dynamical matrix and the Np particles in each of 
these modes with the largest polarization vectors. The 
values of Nm and Np are eventually varied to maximize 
correlation with particle rearrangements. Fig. [2][a) illus- 
trates the locations of the particles identified by the low- 
est Nm = 30 modes and Np — 20 particles for a particular 
configuration. Note that the largest polarization vectors 
are spatially clustered into regions, and that the same 
regions appear in several different modes. Each of the 
Np particles in each of the Nm modes is then assigned a 
value of unity, while the remaining particles are assigned 



a value of zero. We separate this binary map into local- 
ized clusters or "soft spots" as shown in Fig. ^h) . 
Thus, the population of soft spots at strain 7 is repre- 
sented by a binary vector 8(7) = {Si{j) G {0, 1}}, where 
S'i = 1 if particle i is in a soft spot and Si = other- 
wise. In addition, we construct a binary vector for each 
soft spot, indexed by a: Sa = {sa,i{j) G {0,1}}, with 
Sa,i = 1 if particle i is in soft spot a and Sa,i = 0, oth- 
erwise. Thus, 8(7) = J2a^ai^)- Note that to calculate 
the soft spots we used only structural information (the 
particle positions and interactions). As a result, the soft 
spots are structural, not dynamical features. 

We now calculate the correlation of each soft spot, 
80(7) with the next rearrangement at strain 7c, R(7c) = 
{Riilc)}, where i?^ = 1 if particle i has one of the 
largest displacement vectors in the critical mode and 
Ri ~ otherwise. Here, Ua is the number of particles 
in soft spot a. The correlation is [29] 



s„(7)-R(7c) , (l-s,,(7)).(l-R(7e)) 
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The correlation is unity if and R are perfectly cor- 
related, and zero if they are uncorrelated. The rearrange- 
ment shown by the red arrows in Fig. [2] has Cf — 0.64 
with the blue soft spot (the "best" soft spot with the 
highest value of C""", which we label as a = 1). 

The correlation depends on the number of modes, 
Nm, and number of particles per mode, iVp, used to define 
the soft spots. We choose A'^ and TVjJ^ to maximize the 
correlation Cf (Eq. [T|) with the best soft spot, averaged 
over all strains studied. We find N* — 20 particles per 
mode with N'^ = 30 modes, corresponding to roughly 13 
soft spots in a 2500-particle system; however, the results 
are not too sensitive to Np and Nm near their maximum 
values as long as the fraction of particles in soft spots is 
fixed at = 0.1. 

What is the physical significance of the maximal values 
N^ and N*! To understand why the N;^ = 30 lowest 
frequency modes are singled out, we examine the distri- 
bution of polarization vector magnitudes for each mode. 
Each normal mode is composed of N d— dimensional po- 
larization vectors that specify the displacement of each 
particle in the packing. Fig. [3] shows polarization vector 
distributions for (a) the 15 lowest frequency modes and 
(b) 50 intermediate frequency modes. The location of 
these regimes are indicated on a plot of the density of 
states, D{uj) in the inset to Fig.[3l^a). 

Fig. ^ih) shows that for modes in the middle of the 
spectrum, corresponding to extended anomalous modes 
of the type described by Wyart, et al. 0] that con- 
stitute the boson peak the distributions appear to 
be universal with a form given by a modified Gaus- 
sian Orthogonal random matrix ensemble (solid line in 
Fig-EUb)) [la]- While most of the modes in the spectrum 
are well-described by this universal curve, there are clear 
deviations at the low and high frequency ends of the spec- 
trum. At the high frequency end, the localized modes 
differ from the universal curve but play no role in our 
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FIG. 3; (color online) Polarization vector magnitude distri- 
butions for normal modes of a 2500 particle system, (a) The 
15 lowest frequency modes, and (b) 50 "extended anomalous" 
modes from the middle of the spectrum. The dashed line rep- 
resents the distribution for a random matrix ensemble. Inset: 
Density of states D{oj) as a function of frequency showing the 
frequency ranges of the modes in (a) and (b). 



analysis. At the opposite end of the spectrum, Fig. ^a.) 
shows that the lowest-frequency modes (0(10) for a 2500- 
particle packing) also differ significantly from the univer- 



sal curve (11[. This number is consistent with ~ 30 
and a different analysis by Schober and Oligschleger [2lj . 

The quantity N* is the number of particles per mode 
in soft spots. We estimate the size of an individual 
spot independently by analyzing the number of particles 
that change neighbors during "elementary" particle re- 
arrangements. In our quasi-static simulations, "elemen- 
tary" particle rearrangements are defined as those where 
the critical mode is at least 80% correlated with the to- 
tal displacement of all the particles after the packing has 
reached a new mechanically stable state. The results are 
not sensitive to the particular threshold used as long as 
we exclude avalanches, in which one rearrangement trig- 
gers another, and so on. We find that the average number 
of particles that change neighbors during an elementary 
rearrangement is 10. This is of the same order of magni- 
tude as Np ~ 30. Thus, and TV* correspond to the 
number of low-frequency modes that differ significantly 
from the extended anomalous modes, and the size of a 
localized rearrangement, respectively. 

Now that we have identified soft spots, we need to 
show that they are good candidates for structural defects, 
analogous to dislocations in crystalline solids, that con- 
trol flow. The following properties of dislocations ensure 




FIG. 4: (a) Correlation of individual soft spots with the re- 
arrangement field as a function of how much additional 
strain is required to initiate a particle rearrangement (7c — 7) 
for the "best" spot with greatest overlap (solid circles), and 
the second (open squares) and third-ranked (open triangles) 
spots, (b) Correlation of soft spot distributions as a func- 
tion of the difi'erence in strain between the distributions, ^7 
(Eq. [2} . The vertical dashed line indicates the average strain 
between particle rearrangements, showing that soft spot dis- 
tributions are correlated across multiple rearrangements. 



that they control flow: (1) rearrangements tend to occur 
at dislocations, (2) dislocations are long-lived compared 
to the time between rearrangements, and (3) dislocations 
are structurally distinct from the rest of the system. We 
now show that soft spots also possess these qualities. 

(1) Rearrangements occur at soft spots. Each rear- 
rangement is much more strongly correlated with one soft 
spot (the "best" one) than any of the others (Fig. UJ^a)). 
Thus, each rearrangement occurs at one and only one soft 
spot in the population. Moreover, the correlation with 
the best soft spot is high even when the spots are iden- 
tified far in advance of the rearrangement. In Fig. HJa) 
the solid symbols show the correlation between the rear- 
rangement and the best soft spot as a function of the dif- 
ference, 7c— 7, between the strain at which the rearrange- 
ment occurs, 7c, and the strain at which the soft spot was 
identified, 7 < 7c. The correlation decays slowly with 
increasing 7c — 7; the best soft spot calculated shortly 
after a rearrangement still has a strong correlation with 
the next rearrangement. 

(2) The population of soft spots is long-lived compared 
to the interval between rearrangements. To calculate the 
correlation Css between soft spot distributions, we first 
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normalize the soft spot distribution so that it has zero 
mean and unit variance: S — {5*^} for all particles z, 

where ^ = (5, - S,)l ^/Y'W^'W^'W ■ We then 
define 

a. (<57) = — dj S (7) • S (7 + <57) , (2) 

Itot Jo 

where jtot is the total strain studied. Fig. Hljb) shows 
that Css decays slowly compared to the average strain 
between rearrangements, 7. The decay strain is approx- 
imately the product of 7 and the total number of soft 
spots (of order 10), consistent with our observations that 
each rearrangement destroys a soft spot. This would im- 
ply that Css decays at a finite rate even in the thermo- 
dynamic limit where 7^0. 

(3) The population of soft spots is structurally different 
from the remainder of the packing. Compared to the rest 
of the system, soft spots have an average coordination 
number that is 6 ± 1% lower, a bond orientational order 
that is 32 ± 4% lower, and an excess free volume per 
particle that is 103 ± 19% higher. However, we could 
not identify the same population of soft spots by coarse- 
graining these geometric quantities over the area of an 
average spot. Thus, although soft spots are structurally 
different from the rest of the system, the difference is 
sufhciently subtle that one cannot identify them correctly 
using only these local geometric quantities. 

We also calculated the local shear modulus [111, aver- 
aged over the area of a soft spot. This structural quantity 
depends on both the local geometry and the interactions. 
As first noted in [22], a sudden drop in the local shear 
modulus predicts the timing and location of the subse- 
quent particle rearrangement. Thus, if one is sufficiently 
close to the rearrangement, the spatial distribution of the 



coarse-grained shear modulus pinpoints when and where 
the rearrangement will occur. However, it does not pro- 
vide information about other soft spots that are suscep- 
tible but do not rearrange. 

The soft spot analysis, on the other hand, provides 
fundamentally different information about the system. 
Unlike the local shear modulus, the soft spot analysis 
identifies a collection of spots; the next rearrangement 
will occur at one of these spots but the analysis does not 
single out that particular spot a priori. 

So why is it useful to identify a population of soft 
spots? The advantage becomes apparent when one con- 
siders the effects of fluctuations that arise from temper- 
ature or shear. We expect that when fluctuations are 
present, a rearrangement will not necessarily occur in 
the spot with the lowest energy barrier, but could occur 
in any one of the spots with some probability. In that 
case, a statistical description of the soft spot population 
and rearrangements is necessary. 

Our results show that the soft spot population, un- 
like the low-frequency vibrational modes from which it 
is derived or regions of low local shear modulus, is long- 
lived compared to the interval between rearrangements. 
This implies that the population is robust to the fluctu- 
ations that arise in quasistatically-sheared systems. Re- 
cent experiments on thermal colloids show that soft spots 
are also robust at nonzero temperatures [2^. Taken to- 
gether, these results provide strong evidence that the 
spots do indeed constitute the structural defects relevant 
for flow in disordered solids. 

We thank S. R. Nagel and Ke Chen for instructive 
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